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We present preliminary results for the T spectrum on an anisotropic lattice using the improved 0{mv^) NRQCD 
Hamiltonian. We find that accurate results can be obtained in moderate computer times and that they agree 
with earlier results on an isotropic lattice. 



1. The anisotropic lattice 

On an anisotropic lattice an improved action is 
used in which the coupling in the "time" direc- 
tion is tuned so that the temporal lattice spac- 
ing at is much smaller than the spatial lattice 
spacing Os, thus allowing refined measurements 
of mass differences whilst allowing a coarse spa- 
tial lattice. The bare anisotropy x is a parameter 
in the action which determines the renormalized 
anisotropy XR = 0's/o,t- We use the improved 
action given by [ |l|, : 
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Here s, s' run over spatial directions, Ps,s' is a 
1x1 plaquette, Rs's',s is a 2 x 1 plaquette, and 
Us and Ut are tadpole improvement parameters 
which are determined self-consistently. The value 
ofXB. is measured, for example, by computing the 
static quark potential in both the coarse and fine 
directions and we suggest below an alternative 
method to measure xr- We use lattices whose 
parameters were given to us by Alford et al.: (3 = 
1-8, X = 4, Xfl = 3.815(10) ,us = 0.7255 = 
0.9812, where the tadpole coefficients are given 
by the respective mean link values in the Landau 
gauge. Typical lattice sizes used were 8^ x 40 and 
6^ X 60. 
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Whilst we have worked with the values given 
above for Ug and Ut we have also implemented 
a fourier accelerated self-consistent calculation of 
the Landau gauge definition of these parameters. 
The definition used for Landau gauge is to maxi- 
mize 

J2 ^ReTr \u^{x) - -^U^{x)U^{x + fi)] 

with respect to gauge transformations, and define 
Mm = {Uii)Landau- this proccdurc is very time- 
consuming. For example, on 8"^ x 32 to fix the 
gauge to 1 in 2 • 10~^ requires between 50 and 
200 iterations taking, on average, 30 minutes on 
a 533 MHz DEC PC. To obtain better than 1% 
accuracy requires about 50 such measurements. 
We have made preliminary determinations of Us 
and Ut on 8^ x 32 and 4"^ x 16 and find a lattice 
size dependence: 

4^ X 16 : Us = 0.680(2) , ut = 0.9792(3), 
8^ X 32 : Us = 0.707(3) , Ut = 0.9811(2). 

The 8'^ values are close to those of Alford et al. 
and the discrepancy should have little effect on 
our results. 

2. NRQCD to 0{mv^) 

In NRQCD the inversion problem of the 
fermion matrix is an initial value problem. We 
use the evolution equation along the Euclidean 
time-direction defined by: 

G{x,t + l-y) = 

U{n,t) (I- atSH) G{x,t;y) , t>0 
Gix,l;y)^Uin,0) S{x,v) , 



(1) 



Here S{x,y) is the source at the first timeshce, 
t = . We have S{x, y) = 5^'^\x, y) for a single 
quark at the origin but we also evolve a quark 
with a smeared source centred at y. For smearing 
on a length scale I the extended source is defined 

by 



S{x,y-l) = 1 + 
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where D is the covariant derivative on x. Typi- 
cally, 1 = 1, TO = 10 . To construct P-wave me- 
son states we also need single quark propagators 
evolved from a source of the form 

P^{x,y;l)=[l + ^\ D,S^^\x,y). (3) 
The NRQCD Hamiltonian to 0{niv^) is given 

by [|i 
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The fields E and B are derived from the im- 
proved field strength tensor 

where Fjj^u is the standard clover definition. We 
set Ci = 1, i = 0, . . . , 8 . 



In this preliminary study we calculate the low- 
lying S'-wave triplet and singlet states P-wave sin- 
glet states. 

For 5-wave states we use the propagators 

singlet: rf(p,i;y) = 

cos (p • a;)ReTr (Gf t (x, t; y)Gf{x, t; y)) , 

x 

triplet: Tf{p,t;y) = 
J2 cos (p ■ x)ReTT (crGf t {x, t; y)(TGf{x, t; y)) , 

X 

where is evolved with a source defined by eqn. 
1^, where p is an appropriate lattice momentum 
and where the subscripts label different initial 
sources. In practice, we considered only a = 1 
corresponding to the local delta-function source 
and a = 2 corresponding to smearing by / = Og. 
For P-wave states we use the propagator 

^ReTr(Gfi(a;,t;y).AGf(a;,i;y)) , 

X 

where G^ has been evolved from a source defined 
by eqn. ^ with label i denoting the spatial deriva- 
tive. Again, we choose a = 1, 2 defined as for the 
S'-wave states. The averaged propagators are fit- 
ted by a full two- or three-exponential fit using 
the standard SVD technique with full correlation 
matrix to give the various mass differences. The 
value of Of is determined from the l^Pi — l^Si 
mass difference using the experimental value of 
440MeV and then = at/xR- 

In the case of the state the kinetic mass, 
Mfc, is determined by fitting the exponent of the 
leading exponential term to 

Es{p) = -Bs(O) + //2Affe + [p'f/ml . 

The natural unit in the simulation for is 
flf/a^ — l/(atx|j). The bare quark mass enter- 
ing the definition of H, eqn. ^, is chosen so that 
this kinetic mass agrees with the experimental T 
mass. In terms of the computed dimensionless ki- 
netic mass, Mfc, we have Mk = Mk/o,tXR ■ This 
equation yields a method for measuring the renor- 
malized anisotropy xr. Given the dimensionless 
kinetic masses for the these states, and 
respectively, we find 



Xl = iM - Ml)l{Ef 



Table 1 



L nibat 




(liPi - l35i)a* 


a^' MeV 


a-^MeV 


MkX^at 


Mk MeV 


6 1.5 




0.1453(8) 


3028(16) 


794(5) 


50.05(16) 


10413(100) 


8 1.5 




0.1467(40) 


2999(80) 


786(20) 


48.4(1) 


9973(300) 


8 1.7 




0.1584(40) 


2780(70) 


723(20) 


58.5(2) 


11174(250) 


8 1.9 










68.7(2) 


13122(250) 


Table 2 




L 


AEat 


AE MeV 


A£; (isotropic) MeV 


A£:(expt.) MeV 


1^5*1 — 1^5*0 


6 


0.00598(2) 


18.11(7) 


26.1(1) 






8 


0.0060(3) 


18.0(1.0) 






2^5*0 — 1^5*0 


6 


0.190(4) 


575(12) 


610(50) 






8 


0.185(4) 


554(20) 






235'i - l^Si 


6 


0.185(3) 


560(10) 


560(50) 




560 




8 


0.186(4) 


558(20) 






liPi - 13^1 


6 


0.1453(8) 


440 


440 




440 




8 


0.1467(40) 










2iPi - l^Si 


6 


0.37(1) 


1120(33) 






775 




8 


0.38(1) 


1140(33) 









where Ep and Es are the dimensionless expo- 
nents of the leading decay in each case. 

The gauge field configurations were generated 
by using Cabibbo-Marinari and microcanonical 
updates. The simulation was carried out on the 
Hitachi SR2201 parallel computer at the Univer- 
sity of Cambridge High Performance Computing 
Facility. An 8 hour job produced about 2500 in- 
dependent configurations and propagator. 

3. Results and Conclusion 

In Table we show the results on a x T 
lattice for the l^S'i kinetic mass as a function of 
•nih (eqn. ^) and the determination of the lattice 
spacing from l^Pi — l^^i. We have not com- 
puted the spectrum for rrif, = 1.9 but the lattice 
spacings should change only little from the val- 
ues quoted for the other masses. We worked with 
mh — l.ba^^ for which Mk is a bit high but tol- 
erable. Our value for aj^ ^ SOOMeV agrees well 
with the value communicated to us by Alford et 
al. The S and P state mass differences in the 
various cases arc given in Table ^ In the sec- 
ond to last column are the results from our high 
statistics analysis of UKQCD quenched isotropic 



lattice configurations on 16^ x 48 at /3 = 6.0 [ 
^. A combined fit to the l^Si and 1^5*0 data 
gives the triplet-singlet mass difference. The fits 
to both S and P wave data are very good. 

Our preliminary study of NRQCD using an im- 
proved action on an anisotropic lattice has shown 
that the results for heavy quark states agree well 
with similar studies using a standard action on 
isotropic lattices of smaller spatial lattice spac- 
ing, and the accuracy with which parameters can 
be measured with modest amounts of computer 
time is superior. There is no lattice size depen- 
dence in changing from L = 8 to L = 6. 
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